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The method involves preliminary triangulation for 
modelling of boundaries, and division of the 
region into blocks with up to six faces 
conresponding to discontinuities or principal 
faults. Each face is reticulated by interpolation 
from its edges to form an intermediate grid, 
whose nodes are projected from a defined 
direction on to the triangulated surface. Each 
block is stocked with a regular distribution of 
internal nodes by interpolation from those of the 
grids formed on the various faces. The 
positioning of nodes of the elementary cells is re- 
arranged by relaxation in accordance with 
predefined criteria. 
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Description 

[0001 ] La pr6sente invention concerne une m^thode pour g§n6rer un nnalllage 3D dans le but de mod^liser un corps 
quelconque dont la geom6trie est connue et innpos6e a priori tel que par exemple une structure gdologique en tenant 

5 compte fidelement de sa g§ometrie dans toute sa complexite naturelle. Le nnaillage realise par la methode s'accorde 
bien dans le cas d'une structure g^ologlque, avec la position des interfaces entre les couches, avec les failles meme 
de fonne complexe, les biseaux et autres irregularites de structure, en les representant avec leur confonnation reelle. 
[0002] La methode selon {'invention trouve des applications notamment pour la modelisation en geonnecanique, pour 
la caract^risation des reservoirs : couches de sediments, stratigraphie, et aussi pour la simulation des dcoulements 

10 dans les bassins et reservoirs. 

I) ETAT DE lA TECHNIQUE 

[0003] La modelisation d'une structure geologique est rSalisee k partir d'un ensemble de donn^es disponibles ac- 
ts quises soit par des mesures, soit par des analyses et des interpretations de donn^es sismiques ou des diagraphies 
realisees dans des puits. Ces donnees se pr^sentent sous la forme d'un nuage de points entre tesquels, au mieux, il 
convient d'installer un maillage. La formation d'un maillage en 3D joue un role crucial pour faire des simulations nu- 
m^riques dans differentes applications gdologiques. Beaucoup de recherches ont et6 menses sur le sujet principale- 
ment en ingenierie de reservoir en raison de ses applications dans le domaine petrolier. 
20 [0004] On connait de nombreuses techniques de maillage 3D pour des simulations numertques par elements ou 
differences finis et des ensembles de logiciels de simulation cr4es pour leur mise en oeuvre qui sont adaptes k des 
applications et respectent des contraintes detemninees. 

[0005] II est conn u pourcertaines applications d'utiliser des maillages de type tetraedrique non structure, tl est connu 
aussi de faire du maillage hexa^drique de peau ou de surface et de r^partir des points dans le maillage soit par 

25 interpolation transfinie, soit encore par resolution d'equations aux deriv^es partielles. Les surfaces sont generalement 
cr66es k partir de fonctions parametriques et les points du maillage sont obtenus simplement par §chantiilonnage de 
ces fonctions. Une telle methode ne convient pas generalement pour des applications dans les geosciences ou les 
volumes a modeiiser ont des formes le plus souvent tres complexes, difficiles a traduire par des fonctions parametriques 
compte-tenu de la precision qui est recherchee. Aussi prefere t-on approcher les formes complexes par des elements 

30 de surface trianguies. 

[0006] Differentes techniques de maillage appliquees a la modelisation de reservoir, sont decrites par exemple dans : 

Georges P.L. : "Generation Automatique de Maillages", Editions Masson, 1991 ; 

35 . Guerillot D.R. et al: An ineractive 3D Mesh Builder for Fluid Flow Reservoir Simulation; in SPE - Computer Appli- 
cations, Dec. 1993; 

Aziz K.; Reservoir Simulation Grids : Opportunities and problems. SPE-JPT, July 1993; ou 

40 - Heineman Z.E. et al : Gridding Techniques in Reservoir Simulation. IntI Forum on Reservoir Simulation. Alpbach 
Austria, Sept. 12-16, 1988. 

[0007] Cependant, la plupart des techniques de maillage connues, et de 1^ les simulations qui en resultent, ne con- 
viennent bien que pour modeiiser des stmctures geologiques relativement simples. Les irregularites geologiques cou- 

45 rantes dans les zones-reservoir modelisees, telles que les couches ayant de grandes variations rfepaisseur, les failles 
ou les biseaux, sont rarement pris en compte ou alors seulement sous une forme grandement simpllfiee. Le document 
RENARD P ET AL: 'Three-dimensional geometric modeling of a faulted demain: the Soultz Horst example (Alsace, 
France)' COMPUTERS & GEOSCIENCES, NOV. 1994. UK, vol. 20, no. 9, pages 1379-1390, prend en compte les 
irregularites geologiques en prevoyant un traitement specifique des volumes des couches presentantdes irregularites 

50 par construction d'une surface trianguiee qui passe au mieux par les donndes existantes et dont la g6ometrle s'accorde 
avec la configuration geologique. 

[0008] Dans les modeles mailles obtenus dans le cadre des geosciences, les failles apparaissent verticalement; au 
mieux, elles sont inclinees et planes. Or il est important de rendre fidfelement, la geometric en 3D des sediments car 
lis possedent des caracteristiques physiques propres, et aussi la geometrie des failles et autres disco ntinuites meca- 
55 niques. 

[0009] La methode selon T invention permet de generer un maillage 3D respectant la geometrie d'une zone d'un 
milieu complexe dont la geometrie est connue et imposee a priori, tel qu'une structure geologique, dans le but de 
reallser un module representatif de ce milieu. La structure est definie k partir de donnees disponibles acquises soit 
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par des mesures, soit par des analyses et des interpretations de donnees sismiques ou des diagraphies realisees 
dans des puits. Elle comporte una modelisation pr6alable des limites (bords de la zone, discontinuites ou failles par 
exemple) par des surfaces triangul6es, enveloppant au mieux ces linnttes, et un fractionnement de la zone d'6tude 
suivant en blocs ayant au plus six faces, dont les faces correspondent k des discontinuites ou des failles principales 
5 de la zone prealablement triangul6es. 

[0010] Elle est caract6rls§e en ce que'elle comporte : 

la subdivision de chaque bloc en un ensemble de cellules 6lementaires necessaires pour modeliser la zone par 
Elements finis, par une reticulation en 3D de chacun des dlff6rents blocs obtenus, subdivision qui comprend : 

10 

a) le maillage de chacune des faces du bloc par interpolation k partir de ses bords, de fapon a cr§er une grille 
intenm6diaire, suivi d'une projection sur la surface triangul6e des noeuds de cette grille interm§diaire suivant 
une direction definie (perpendiculairement h elle-meme par exemple); 

15 b) le peuplement en noeuds de grille de chaque bloc dans son 6paisseur, par une m§thode d'interpolation k 

partir des noeuds des differentes grilles form§es sur les differentes faces triangulees de chaque bloc, et 

c) un rearrangement dans le positionnement des noeuds des cellules 6l6mentaires cr66es, obtenu (de fagon 
iterative, gen^ralement) par relaxation dans le positionnement des noeuds, en fonction de criteres predefinis 
20 tels que la r^gularite, le choix de la tailie relative des differentes cellules et le mode de progression de leurs 

tallies, la direction de projection etc. 

[0011] Chaque bloc etant ainsi reticule en 3D, on peut proceder a un rdassemblage ou reempilement de tous les 
blocs. 

25 [001 2] On peut proceder successivement a une premiere operation de relaxation dans le positionnement des noeuds 
sur chacune des differentes faces triangulees, le peuplement en noeuds Intemes etant effectue k partir des grilles 
relaxees, et ensuite k une deuxifeme operation de relaxation dans le positionnement de ces noeuds Internes. 
[0013] De preference, la premiere operation de relaxation est effectuee par des phases iteratives de relaxation- 
projection. 

30 [0014] La methode selon Tinvention pennet d'obtenir un modfele qui respecte trfes pr6cisement la geometric souvent 
complexe de structures geologiques telles que des champs petrolif feres, ce qui a pour effet que les simulations qui 
peuvent etre faites k partir de ce modele rendent plus fidetement compte des phenomfenes. La methode utilisee penmet 
en outre une economietres importante de temps de calcul pour realiser les operations de maillage 3D habituellement 
tres longues avec les methodes existantes. 

35 [0015] D'autres caracteristiques et avantages de la methode selon {'invention, apparatlront k la lecture de la des- 
cription ci-apres de modes de realisation decrits ^titre d'exemples non limitatifs, en se ref6rant aux dessins annexes ou : 

les Fig. 1 et 2 montrent des exemples de relaxation -projection sur un maillage surfacique 4X4; 

40 - la Fig.3 montre un cube deforme avec des faces triangulees; 

la Fig.4 montre I'enveloppe externe d'un cube mailie; 

les Fig.5 et 6 montrent les surfaces internes d'un cube mailie respectivement avant et apres relaxation; 

45 

la Fig. 7 montre des courbes relevees sur cinq coupes d'une structure geologique de tres grandes dimensions, sur 
laquelle la methode selon Tinvention a ete testee, qui constituent une partie des donnees initiales dont on dispose, 
sur cette structure; 

50 - la Fig. 8 montre d'autres donnees initiales sur la meme structure, en I'occurrence, des cartes d'isobathes et les 
emplacements de 4 puits; 

la Fig. 9 montre une modelisation surfacique ( base de surfaces triangulees) de la structure geologique etudiee; 
55 - la Fig. 1 0 montre un exemple de decoupage de la structure en blocs; 

la Fig. 11 montre le resultat d'un maillage bloc par bloc de I'enveloppe exteme de la structure geologique; 
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la Fig. 12 montre !e r^sultat d'un maillage bloc par bloc des surfaces internes de la structure geologique; 

la Fig. 13 illustre k titre d'exemple comment on prend en compte un biseau dans le maillage de la structure geo- 
logique; 

5 

la Fig.1 4 illustre k titre d'exemple comment on prend en compte une faille dans le maillage de la structure geolo- 
gique; 

la Fig.15 illustre la m§thode d'interpolation 2D de Coons; 

10 

la Fig.1 6 montre les elements pris en compte dans la methode d'interpolation 3D de Coons; et 

la Fig. 17 illustre le mode de relaxation-projection des noeuds de maille. 

15 [0016] La structure geologique de grandes dimensions illustree par les Fig. 7- 10. con-espond k un dome anticlinal 
comportant six interfaces entre des sediments. 

[0017] La realisation d'un maillage en 3D par elements finis d'une zone a partlr d'un ensemble de donnees initiales 
connues pour la zone : interfaces entre diff§rentes couches ou horizons, failles principales, carte en courbes de niveau 
avec les emplacements de puits passe par une succession d'etapes qui seront decrites ci-apres. 

20 

!) Modeiisation des limites ou bords de la zone etudiee et de ses principales discontinuitds. 

[0018] On realise une modelisation des surfaces : bords, interfaces, failles en uttlisant une logiciel specialise connu 
tel que GOCAD (marque depos^e) sous la fomne d'une juxtaposition continue d'eiements de surface triangulaires. 

25 Dans un premier temps on fomne une ebauche qui s'adapte au mieux aux donnees initiales. De maniere iterative, on 
deplace les sommets de la surface en utilisant un logiciel d'interpolation et on densifie localement les elements de 
surface triangulaires, de faQon k plaquer au mieux la surface sur le semis de points (donnees d'entree). Les elements 
de surface ainsi crees ne se raccordent generalement pas comme lis le devraient : les interfaces ne sont pas limitees 
aux intersections et au bord du domaine d'etude et elles peuvent se crolser en debordant. On calcule alors les lignes 

30 d* interception entre surfaces et on limrte celles-ci aux courbes d'intersection. La Fig 7 montre le r6suttat final de cette 
modelisation des limites de la structure et Ton peut remarquer que la plupart des triangles des surfaces trianguiees du 
biseau sont localises dans les regions de forte courbure. 

[0019] A ce stade, 11 n'estpas encore possible de produire un maillage 3D tenant compte des contraintes geologiques 
et numeriques. Pour y parvenir, on procede a un fraction nement de la zone. 

35 

II) Fractionnement de la zone etudiee 

[0020] On decoupe d'abord la zone en macro-blocs ayant au plus six faces. II s'agit generalement de blocs de forme 
hexaedrique (a six faces), mais le nombre de faces peut etre inferieur a six et descendre jusqu'a deux. Dans ce cas, 

40 pour les besoins de la methode, le bloc est quand meme considerere comme un block hexaedrique mais ayant des 
faces degenerees (reduites k un point ou k une courbe). Le mode de decoupage est generalement suggere par la 
geometrie des horizons, des principales failles et leur interaction. Pour constituer ces macro-blocs 3D, on a besoin 
parfois d'introduire des frontieres artificielles. C'est necessaire parexemple quand un bloc a plus de six faces naturelles 
et ne satisfait done pas k la rfegle choisie d'etre au plus hexahedrique. Dans un tel cas, on reunit deux faces en une 

45 face complexe pour satisfaire le critere, ou bien on introdurt des surfaces de coupure artificielles. 

[0021] Un exemple de decoupage de la structure geologique etudiee en macro-blocs est montre k la Fig. 1 0, ou Ton 
a extrapoie la gaille jusqu'^ la frontifere de la structure et verticalement projete I'intersection du biseau sur les autres 
surfaces. 

50 III) Maillage 

[0022] On realise une reticulation ou maillage en 3D separement de chacun de ces macro-blocs de fa^on k le sub- 
diviser en blocs eiementaires qui vont fonner les cellules pour la modelisation par elements finis recherchees. Cette 
operation est realisee en deux temps. 

55 

1) Maillage des faces 

[0023] a) on cree tout d'abord une grille sur chacune des cinq ou six faces de chaque macro-bloc. Pour cela, on 



4 



EP0 008 0 1364 Fhttp://wv w .ge tt hepat e ntx om/Log in.dog/$ mdean / F etch/ Defa ult .do^^^ 



Page^of 25 



EP 0 801 364 B1 

forme une grille intemnediaire par interpolation a parti r des quatre bords de la face, en utilisant la methode rfinterpolation 
transfinie de Coons dont le principe est rappele ci-apres. 

[0024] Etant donn6 quatre courbes param§triques f1(u), f2(u), g1{v), g2(v) (0 2 u, v 2 1) d^finissant quatre bords 
jointifs (Fig. 15), on introduit une fomaule pour calculer la surface la plus tendue qui passe par les quatre bords (qui 
5 interpole ces bords) : 



5(«,v) = (l-«)gl(vO + u52(v) + (l-v)/l(«) + v/2(u)-[(l-«3(l-v)P(0,0) ^ 
10 +(1 -iOvP(0.1)H-a(l - v)P(1.0)+uvP(Kl)] ^ 

[0025] Les sonnnnets P(i, j) sont les quatre coins de la surface. 

[0026] On peut alors obtenir un maillage regulier NxM de cette surface par simple §chantiHonnage en prenant pour 
^5 sommets les points S(i, j) correspondant aux ui = i/(N-1 ), vj i=0 & N-1 , j=0 ^ M-1 . La formule 1 restreinte au 

points du maillage devlent : 

-[(I - TriqXl - TTTi )Pm) + (1 - Tjfrj PfO, I) (2) 

25 

[0027] Dans le cas illustre, les courtaes f1 , f2, g1 , g2, sont obtenues & partir du logiciel GOCAD pr6crt6 et sont done 
des iignes brisees. Pour obtenir un maillage regulier MxN, on echantillonne regulierement (en abscisse curviligne ou 
longueur de segment) les couriDes f 1 et f2 en N points, et les courbes g1 , g2 en M points, puis on applique la fonmule 
30 (2). On peut tout aussi bien obtenir un maillage k progression geometrique sur chaque direction, simplement en im- 
posant cette contrainte au niveau du reechantillonnage des courbes. 

[0028] Ce maillage intemnediaire sauf cas exceptionnel ne convient pas car il n'y a pas de raison qu'il coincide k ce 
stade avec la surface qui a pr6a!ablement §t6 triangul6e par le logiciel de surfa9age. L'ajustement necessaire consiste 
k projeter les noeuds de la grille Intemnediaire sur la surface triangulee orthogonalement a celle-ci, ou suivant d'autres 

35 directions , de fa9on k obtenir une nouvelle grille avec des noeuds places sur la surface triangulee. Quand les quatre 
courtes frontieres sont regulierement echantillonnees, la methode dMnterpolation de Coons presente I'avantage de 
produire une grille reguliere en ce sens que sur chaque couri^e Isoparametrique tous les segments echantillonnes ont 
la meme longueur. Plus generalement la methode de Coons conserve les rapporfs de taille entre elements consecutifs 
dans une direction donnees. Ceci penmet d'avolr par exemple des maillages regu tiers ou respectant cerfaines pro- 

40 gressions dans la taille des mailles. Cependant la grille obtenue par cette methode connue perd sa regularite (ou le 
critere recherche) une fois qu'elle est projetee sur la surface triangulee. 

2) Rearrangement des noeuds de face 

45 [0029] Pour retro uver le critere desire dans la distribution des noeuds de la grille projetee, une regularite satisf aisante 
par exemple, il faut deplacer les noeuds du maillage sur la surface triangulee. Pource faire, on procede k une relaxation- 
projection de la grille par ajustements success ifs comme indique ci-apres (dans le cas ou le critere est la regularite), 
en relation avec les Fig. 1 , 2. 

[0030] Etant donne un bloc mailie, le but cette relaxation est de rendre son maillage le plus regulier possible au sens 
50 oii n'importe quelle courbe isoparametrique du maillage soit reguliferement echantillonnee. Pour ce faire. on deplace 
de fa^on iterative tout point du maillage de maniere k le positionner au mieux a egale distance de ses deux voisins 
surchacune des trois directions (u,v,w) de voisinage. Soit P le point qu'on cherche a deplacer. Pour chaque direction 
i, (i = u, v.w) (Fig.1 7), on trouve le point Pi qui se situe k egale distance des deux voisins de P sur la courbe de direction 
i. (i = u, v,w). Le point P est alors rem place par le barycentre des trois points Pi (i = u, v, w). 
55 [0031] Aprds avoir deplace tous les points du maillage de la sorte, ceux-ci ne sont plus sur la surface triangulee. On 
les projette alors k nouveau sur la surface triangulee et on repete le processus de relaxation-projection jusqu'^ atteindre 
la regularite souhaitee. 

[0032] Ceci est illustre sur I'exemple des Fig. 1 , 2, correspondant k une face du cube deforme de la Fig, 3 qui est 
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un exemple synth^tique penmettant de bien illustrer les diff^rentes etapes de la methode. Cette surface a ete obtenue 
par simple deformation d'une face carree en preservant ses bords ( le bord de la surface reste un carre). Sur ces 
figures, on peut voirtrois maillages 4X4. Le premier est un maillage plat issu de I'lnterpolation de Coons appliquee 
aux quatre bords (formant un carre). Le deuxieme est le maillage fin dont les noeuds sont obtenus par projection de 
5 ceux du premier maillage. On constate que ce maillage a perdu de sa regularite. Le troisieme est le maillage regulier 
trac6 en gras et dont les noeuds sont situ^s sur la surface. II a et6 obtenu k partir du deuxifeme maillage aprfes une 
serie de projections et relaxations. La Fig. 4 montre le resultat de ces operations appliquees k Tensemble des faces 
du cube defomne. 

10 3) Peuplement en noeuds de grille de chaque bloc dans son epaisseur : 

[0033] Pourprolonger en profondeur les grilles form6es sur les faces des blocs (Fig. 5), on utilise une generalisation 
au volume de la methode d'interpolatlon transfinie pr^cedente. 

[0034] Etant donne 12 courbes parametriques f1 (u), f2(u ), f3(u). f4(u), g1 (v). g2(v), g3(v), g4(v), h1(w), h2(w). h3 
15 (w), h4(w) (0 2 u, V, w2 1) definissant les 12 bords d'un bloc tel que montr6 dans la Fig.16, une generalisation de la 
fomnule de Coons fournit une parametrisation volumique a Tinterieur du bloc. La fomnule de Coons dans le cas volu- 
mique est la suivante : 



on 

V(u, V, vv) = (I - v)( I - vv )/!(«) + (l - v)vv/2(u) + vw/3(u)+ v(l - w) f4(u) 
-(i -u)(l - vw)5l(v)+ ( I - u)wg2{\:) + i/W'g3(v)4- /f( I - w}84(v) 
-K 1 - «X I - vO/i if vv) + ( I - u)vh2{\v) + uvh3(w) - w( 1 - v)h4(v) 

25 +C(«.V-. IVJ 



30 



35 



(3) 



avec 



c( uy , w)=-2[(^ - u}{^ -v)(i - w) P(o,o,c)+(i - t/)(i -v) ivP(o,o,i ) 

+(1-u)v(1-w)P(0,1,0)+(1-u)v»vp{0,1.1)+Lf(1-v)(1-iv)P(1,0,0) 
+u(1 -v)ivP(1 .0,1 )+iAf(1 -iv)P(1 ,1 ,0)+yvivP(1 .1 .1 )] 



[0035] Les sommets P(i, j, k) (I, j, k= 0 ou 1 sur la Fig 1 6) sont les 8 coins du bloc. 

[0036] En combinant les equations (1) et (3) on aboutit k une expression de la fonnule de Coons dans le cas volu- 
mique en fonction des 6 surfaces de bord du bloc, S1(u. w), S2(u.w), S3(v.w). S4(v.w). S5(u.v), S6(u,v). Mais id, les 
40 surfaces sont aussi supposees etre obtenues par interpolation de Coons k partir de leurs courbes de bord. La formule 
est la suivante : 



45 V(M, v,u;) = - ((1 - v)51(u, + vS2{u, w) + ( 1 - u)S3(v^\v) + «54(v'. vw) 

4 

50 

[0037] Un echantillonnage quelconque de V(u,v,w) fournit alors un maillage EF du bloc. 

[0038] Cette interpolation en profondeur donnerait une distribution regulidre de noeuds si elle etait appliquee k partir 
des noeuds des grilles fomies sur les differentes faces du bloc. Ce n'est pas le cas et la distribution obtenue est 
generalement irreguiiere car la methode d'interpolatlon en profondeur est appllqu6e aux noeuds reportes par relaxation 
55 sur les surfaces trianguiees, comme decrit dans retape 2 pr6cedente. 
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4) Rearrangement des noeuds de bloc 

[0039] Aussi, pour retrouver una r^gularitd suffisante (ou respecter un autre critere choist), la m^thode selon ('inven- 
tion comporte 6galement une relaxation de la grille 3D de la manlere indiquee precedemment. Les Fig .5, 6 montrent 
5 le maillage interne du cube deforme avant et apres relaxation. Le nnailtage de la Fig.6, on le note, possede une regularity 
bien nneilleure que le maillage de la Flg.5. 

[0040] Quand chacun des blocs precedents de la structure etudiee a ete convenablement maille en volume, on 
reassemble tous ces macro-blocs apres leur reticulation 3D. 

[0041] La m^thode qui vient d'etre decrite peut §tre utilis^e pour mailler des heterog§n41t§s diverses telles qu'un 
10 chenal ou encore pour mailler radialement un puits devi§. 

[0042] On a decrit des modes de realisation de la methode appliquee au maillage d'une structure geologique. On 
ne sortlralt pas du cadre de I'invention, cependant, en utilisant la methode pour le maillage de tout corps ou structure 
dont la geometric est connue et imposee a priori. 

15 

Revendications 

1 . Methode pour gSnerer un maillage 3D dans une structure dont la g^ometrie est connue et imposee a priori telle 
que par exemple une structure geologique definie a partir de donnees disponlbles acquises par des mesures, par 

20 des analyses et/ou des interpretations de donnees sismiques ou des diagraphtes realisdes dans des puits au 

travers de la structure, dans le but de r§aliser un module reprdsentatif de cette structure, comportant une mode- 
lisation prealable des limites par des surfaces triangulees, enveloppant au mieux ces limites, et un fraction nement 
de la zone en blocs ayant au plus six faces, dont les faces correspondent a des disco ntlnuites ou des failles 
princlpales de la zone pr^alablement triangulees, la methode etant caracterisee en ce qu'elle comporte : 

25 

la subdivision de chaque bloc en un ensemble de cellules elementaires n^cessaires pour mod^liser la zone 
par §l§ments finis, par une reticulation ou maillage en 3D de chacun des differents blocs obtenus, subdivision 
qui comprend : 

30 a) le maillage de chacune des faces du bloc par interpolation k partir de ses bords, de fa^on a cr§er une 

grille intermediaire, suivi d'une projection sur la surface triangulee des noeuds de cette grille intermedtaire 
suivant une direction definie, de fa^on k g6nerer des grilles de face; 

b) le peuplement de chaque bloc en noeuds internes par interpolation a partir des noeuds des differentes 
35 grilles fomnees sur les differentes faces triangulees de chaque bloc, de fagon k obtenir une distribution 

reguliere des noeuds dans le volume de chaque bloc; et 

c) un rearrangement dans le positionnement des noeuds des cellules elementaires creees, obtenu par 
relaxation dans le positionnement des noeuds. en fonction de criteres predef Inis. 

40 

2. Methode selon la revendication 1 , caracterisee en ce que I'on precede successivement a une premiere operation 
de relaxation dans le positionnement des noeuds sur chacune des differentes faces triangulees, le peuplement 
en noeuds internes etant effectue a partir des grilles relaxees, et ensuite une deuxieme operation de relaxation 
dans le positionnement de ces noeuds internes. 

45 

3. Methode selon la revendication precedente, caracterisee en ce que Ton realise la premiere operation de relaxation 
par des phases iteratives de relaxation-projection, 

4. Methode selon Tune des revendications 1 ou 2 ou 3, caracterisee en ce I'on choisit comme direction de projection, 
50 une direction perpendiculaire k la grille intemnediaire. 

5. Methode selon Tune des revendications precedentes, caracterisee en ce que Ton utilise un mode d'interpolation 
de type transfini. 

55 6. Methode selon I'une des revendications precedentes, caracterisee en ce que I'on rearrange te positionnement 
des noeuds des cellules elementaires creees de fa9on k obtenir des maillages reguliers. 

7. Methode selon I'une des revendications 1 ^5, caracterisee en ceque Ton rearrange le positionnement des noeuds 
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des cellules 6l6mentaires cr66es de fa^on k obtenir un mode d6tennin6 de progression de la taille des mailles. 

8. M^thode selon Tune des revendications pr6c§dentes, caracterisee en ce que I'on proc^de k un r^assennblage 
de tous les blocs. 

5 

Patentanspriiche 

1. Verfahren zum Erzeugen eines 3D-Netzwerkes in einer Struktur, deren Geometrie bekannt und a priori vorge- 
10 schrieben ist, wie zum Beispiel eine geologische Struktur, die aus verfugbaren Daten definiert ist, welche durch 

Messungen, durch Analysen und/oder Interpretationen selsmischer Daten oder Diagraphien aufgenommen sind, 
welche in Bohrlochern entlang der Struktur realisiert sind, mit dem Ziel, ein fiir diese Struktur reprasentatives 
Modell zu verwirklichen, umfassend eIne vorherige Modellierung der Grenzen durch aus Dreiecken zusammen- 
gesetzte Oberflachen, welche diese Grenzen am besten umhullen, und eine Fraktionierung der Zone In Blocke 
15 mit hochstens sechs Rachen, deren Flachen Hauptdiskontinuitaten oder Hauptverwerfungen der zuvor aus Drei- 

ecken zusammengesetzten Zone entsprechen, wobei das Verfahren dadurch gekennzeichnet ist, dass es um- 
fasst: 

die Unterteilung jedes Blockes in eine Zusammensteilung von Elementarzellen, die notwendig sind, um die 
20 Zone durch Endelemente durch eine Vernetzung oder 3D-Netzwerk von jedem der verschiedenen erhaltenen 

Blocke zu modeliieren, welche Unterteilung umfasst: 

a) Das Netzwerk von jeder der Flachen des Blocks durch Interpolation ausgehend von seinen Grenzen 
derart, dass ein Zwischengitter erzeugt wird, gefolgt von einer Projektion der Knoten dieses Zwischengit- 

^5 tens auf die aus Dreiecken zusammengesetzten Flache gema(5 einer definierten Richtung, derart, dass 

Flachengitter erzeugt werden; 

b) die Population jedes Blockes mit inneren Knoten durch Interpolation aus Knoten von verschiedenen 
gebildeten Gittern auf die verschiedenen aus Dreiecken zusammengesetzten Flachen jedes Blocks, der- 

30 art, dass eine regelmassige Verteilung der Knoten in dem Volumen jedes Blockes erhalten wird; und 

c) Einer Wiederanordnung in derPositionierung der Knoten dererzeugten Elementarzellen, erhalten durch 
Entspannung in der Pos'rtionierung der Knoten in Abhangigkeit vorgegebener Kriterien. 

35 2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, dass man aufeinanderfolgend zu einem ersten Entspan- 
nungsvorgang in derPositionierung der Knoten auf jeder der verschiedenen aus Dreiecken zusammengesetzten 
Flachen ubergeht. wobei die Population mit inneren Knoten ausgehend von den entspannten Gittern durchgefCihrt 
wird und man anschlieBend zu einem zweiten Entspannungsvorgang in der Positionierung dieser inneren Knoten 
ubergeht. 

40 

3. Verfahren nach dem vorhergehenden Anspruch, dadurch gekennzeichnet, dass man den ersten Entspannungs- 
vorgang durch iterative Relaxations- Projektions-Phasen durchfuhrt. 

4. Verfahren nach einem der Anspruche 1 oder 2 oder 3, dadurch gekennzeichnet, dass man als Projektions- 
^5 Richtung eine Richtung wahit, die senkrecht zum Zwischengitter liegt. 

5. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, dass man einen Interpolati- 
onsmodus vom transfiniten Typ verwendet. 

50 6. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, dass man die Positionierung 
der Knoten dererzeugten Elementarzellen derart wiederanordnet, dass regelmaBige Netzwerke ertialten werden. 

7. Verfahren nach einem der Anspruche 1 bis 5, dadurch gekennzeichnet, dass man die Positionierung der Knoten 
der erzeugten Elementarzellen derart wiederanordnet. dass ein vorgegebener Progressionsmodus des Umfangs 

55 der Maschen erhalten wird. 

8. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, dass man zu einem emeuten 
Zusammensetzen aller Blocke ubergeht. 
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Claims 

1. Method for generating a 3D mesh in a structure whose geometry is known and imposed a prion, such as for 
example a geological structure defined from available data acquired by measurements, by analyses and/or inter- 
5 pretations of seismic data or logs produced in bore holes through the structure, with the object of producing a 

representative model of this structure, comprising prior modelling of the limits by triangulated surfaces, enveloping 
these limits as well as possible, and splitting of the zone into blocks having at most six faces, the faces of which 
correspond to previously triangulated major faults or discontinuities in the zone, the method being characterised 
in that tt comprises: 

10 

subdivision of each block into a set of elementary cells necessary to model the zone by finite elements, by 
creating a 3D lattice or mesh from each of the different blocks obtained, said subdivision comprising: 

a) creating a mesh from each of the faces of the block by interpolation from its edges so as to create an 
'5 intemnediate grid, followed by projection onto the triangulated surface of the nodes of this intermediate 

grid in a defined direction so as to generate face grids; 

b) population of each block with intemal nodes by interpolation from the nodes of the different grids f omried 
on the different triangulated faces of each block so as to obtain a unifomn distribution of the nodes in the 

20 volume of each block; and 

c) rearrangement in the positioning of the nodes of the elementary cells created obtained by relaxation in 
the positioning of the nodes as a function of predefined criteria. 

25 2. Method according to claim 1 , characterised by successively a first operation for relaxation in the positioning of 
the nodes on each of the different triangulated faces, the population with intemal nodes being effected from the 
relaxed grids, and then a second operation for relaxation in the positioning of these internal nodes. 

3. Method according to the preceding claim, characterised in that the first relaxation operation is carried out by 
30 iterative phases of relaxation and projection. 

4. Method according to one of claims 1 or 2 or 3, characterised in that a direction perpendicular to the intemriediate 
grid is chosen as the direction for projection. 

35 5. Method according to one of the preceding claims, characterised in that a mode of interpolation of transf inite type 
is used. 

6. Method according to one of the preceding claims, characterised in that the positioning of the nodes of the ele- 
mentary cells created is rearranged so as to obtain unifomn meshes. 



40 



45 



7. Method according to one of claims 1 to 5, characterised in fAiaf the positioning of the nodes of the elementary 
cells created is rearranged so as to obtain a given mode of progression for the size of the mesh. 

8. Method according to one of the preceding claims, characterised in thatoW the blocks are reassembled. 
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Fig. 1 
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Fig. 2 
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Fig. 3 
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Fig. 4 
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Fig. 5 
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Fig. 6 



15 



EP000801364 Fhtt p://wvw.qetthepatentxom/Lo q in.do g/$ mdean/Fetch/Defau(t.dog/EP000801364xpc?fro m Cach 



Page 16 of 25 



EP 0 801 364 B1 




UL 



16 



EP00 0801364 Fhttp:/AAAAW.getthepatentxoiTi/Lo g in.do q/S m d e an/ Fetch/Default . dog/EP00080136 



Page 17 of 25 



EP 0 801 364 B1 




17 



^ P000801364 [http://wvw.Qetthepatentxom/Lo g in.dog/$rTidean/Fetch/Defauit.dog/EP000 80 1364xpc?fromCa 



EP 0 801 364 B1 



Page 18 of 25 




u. 



18 



EP000801364 fhttp://ww w.getthe palent.com/Login.dog/$ mdean/Fetch /D efault .dog/EP000801364xpc? fromCache=1part=m aintoolbar=bottom Page 19 of 25 



EP 0 801 364 B1 




O 

Li. 



19 



E P0Q0S013 64 [http://w ww. g etthe p atentxom/Log i n.dog/$mdean/Fetch/Default .dog / E P00 0801364 



Pa g e 20 of 25 




20 



.^P-QOPSOI 364 rhttp://www.getthepat e nt.com/Logi n.d og/$ mdean/Fetch/Defauitd og/E PQ0080 1 364 .^ 



EP 0 801 364 B1 



Page 21 o f 25 




21 



cP000801364 [h tt p:/AftWw.getthepatent.c om/Login .d o g/$mdean/ Fetch/Default.do g/EP000 801 3 ^ 



Page 22jpf 25 



EP 0 801 364 B1 




CO 
U_ 



22 



EP0008 01364 f http://vww.getthepate nt. c om/Login.dog /$mdean/Fetc h /Default.dog/EP000801364xpc?from 



EP 0 801 364 B1 



Page 23 of 25 




23 



pP 00Q8Q1364 r http://wvw.g ett hepatent.com /Log in.doq/Smdean/F8tch/Default .d og/E P000 80 1364.cpc?fro m 



Page 24 of 25 



EP 0 801 364 B1 




1/3 

d) 

U. 



24 



EP0008 01364[http://vww.getthepatent .com/ LoginJog/$mdean/Fetch/Default.do9/EP000B 01364xpc ?fro Pa3s25. of 25 



EP 0 801 364 B1 



FIG.16 




FIG.t7 




25 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the appHcant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



